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Abstract 

We present a measurement of the cross section for production of two or more jets as 
a function of dijet mass, based on an integrated luminosity of 86 pb -1 collected with 
the Collider Detector at Fermilab. Our dijet mass spectrum is described within errors 
by next-to-leading order QCD predictions using CTEQ4HJ parton distributions, and 
is in good agreement with a similar measurement from the D0 experiment. 
PACS numbers: 13.85.Rm, 12.38.Qk, 
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Hard collisions between protons and antiprotons predominantly produce dijet 
events, which are events containing at least two high energy jets. A measurement 
of the dijet mass differential cross section provides a fundamental test of Quantum 
Chromodynamics (QCD) and a constraint on the parton distributions of the proton. 
We previously reported measurements of the inclusive jet transverse energy (Et) 
spectrum U and the cross section for events with large total Ex @. Both measure- 
ments indicated an excess of events at high E T compared to the predictions of QCD. 
This letter presents our most recent measurement of the dijet mass spectrum and 
compares it with the predictions of next-to-leading order QCD and the measurement 
of D0 |4|]. This measurement, with an integrated luminosity of 86 pb _1 , is signifi- 
cantly more sensitive to events at high dijet mass than our previous measurements 
of the dijet mass spectrum |J with integrated luminosities of 4.2 pb _1 and 26 nb _1 . 
We recently used this data sample combined with 20 pb -1 of older data to measure 
dijet angular distributions |J and to search the dijet mass spectrum for new particles 
decaying to dijets 0. 

A detailed description of the Collider Detector at Fermilab (CDF) can be found 
elsewhere ||. We use a coordinate system with the z axis along the proton beam, 
transverse coordinate perpendicular to the beam, azimuthal angle 0, polar angle 9, 
and pseudorapidity rj = — lntan(0/2). Jets are reconstructed as localized energy de- 
positions in the CDF calorimeters, which are arranged in a projective tower geometry. 
The jet energy, E, is defined as the scalar sum of the calorimeter tower energies inside 
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a cone of radius R = y(A^) 2 + (A0) 2 = 0.7, centered on the jet direction. Jets that 
share towers are combined if the total Et of the shared towers is greater than 75% 
of the Et of either jet; otherwise the towers are assigned to the nearest jet. The jet 
momentum, P, is the vector sum: P = J2 Qui, with being the unit vector pointing 
from the interaction point to the energy deposition Ei inside the cone. The quantities 
E and P are corrected for calorimeter non-linearities, energy lost in uninstrumented 
regions of the detector, and energy gained from the underlying event and additional 
pp interactions. We do not correct for energy lost outside the clustering cone, since a 
similar loss is present in the O(a^) QCD calculation in which an extra gluon can be 
radiated outside the jet clustering cone. The jet energy corrections increase the mea- 
sured jet energies on average by 20% (16%) for 100 GeV (400 GeV) jets. Full details 
of jet reconstruction and jet energy corrections at CDF can be found elsewhere 0. 

We define the dijet system as the two jets with the highest transverse momentum 
in an event (leading jets) and define the dijet mass as M = J (E\ + E2) 2 — (A + A) 2 - 
Our data sample was obtained using four triggers that required at least one jet with 
uncorrected cluster transverse energies of 20, 50, 70 and 100 GeV, respectively. After 
correcting the jet energies these trigger samples were used to measure the dijet mass 
spectrum above 180, 217, 292, and 388 GeV/c 2 , respectively, where the trigger effi- 
ciencies were greater than 97%. The four data samples corresponded to integrated 
luminosities of 0.091, 2.2, 11, and 86 pb -1 respectively. We selected events with two 
or more jets and required that the two leading jets have pseudorapidities of \r]\\ < 2 
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and 1 7/2 1 < 2 and satisfy | cos 6 1 * | = | tanh[(r/i — 772) /2] | < 2/3, where 9* is the scat- 
tering angle in the dijet center-of-mass frame. The cos 6* requirement ensures full 
acceptance as a function of the dijet mass. The z position of the event vertex was 
required to be within 60 cm of the center of the detector; this cut removed 6% of the 
events. Backgrounds from cosmic rays, beam halo, and detector noise were removed 
by requiring Et / y/J2 Et < 6 GeV 1 / 2 and E < 2 TeV, where Et is the missing 
transverse energy | 1QR , Et is the total transverse energy (scalar sum), and J2E is 



the total energy in the event. These cuts selected 60,998 events. 



The dijet mass resolution was determined using the PYTHIA |TT| Monte Carlo 
program and a CDF detector simulation. The true jet is defined from the true Et 
of particles emanating from the hard scattering, using the same jet algorithm as de- 
scribed above, but applied to towers of true Et- The true Et of a tower is the Et 
of the generated particles that enter the tower. The simulated jet uses the E T of 
simulated calorimeter towers and the jet energy corrections for the CDF detector 
simulation. The Et of the simulated jets is corrected to equal the Et of the corre- 
sponding true jet on average. The dijet mass resolution function, p(M,m), is then 
defined as the distribution of simulated dijet masses, M, for each value of true dijet 
mass, m. The dijet mass resolution was determined for six values of m between 50 
and 1000 GeV/c 2 and then a single smooth parameterization was used to interpolate 
between these values. The dijet mass resolution is approximately 10% for dijet masses 
above 150 GeV/c 2 . 

The steeply falling dijet mass spectrum is distorted by the dijet mass resolution. 



We correct for this distortion with an unsmearing procedure. Define the smeared 
spectrum, S(M), as the convolution of the true spectrum, T(m), and the dijet 
mass resolution: S(M) = f T(m)p(M,m)dm. We parameterize the true dijet mass 
spectrum with T(m) = A(l — m/^/s + Cm 2 / ' s) N /m p where y/s = 1800 GeV. We 
then fit the smeared spectrum to our data to find the value of the four parameters 
A = 6.67 x 10 17 pb/(GeV/c 2 ), C = 2.95, N = -6.98, and P = 6.70. The fit has a x 2 
of 20.5 for 14 degrees of freedom. The unsmearing correction factors, Ki, are then 
defined as the ratio of the smeared to true spectrum, Ki = J i S(M)dM/ f i T(m)dm, 
where the integration is over mass bin i. The value of Ki smoothly decreases from 
1.07 at M = 188 GeV/c 2 , to 1.03 at M = 540 GeV/c 2 , and then smoothly increases 
to 1.12 at M = 968 GeV/c 2 . The corrected cross section as a function of dijet mass 
is given by 

da/dM = m/(Ki C a AM), (1) 

where for each mass bin i, rii is the number of events, C is the integrated luminosity, 
6i is the efficiency of the trigger and ^-vertex selections, and AM is the width of the 
mass bin. 

In Table I we list 12 independent sources of systematic uncertainty in the dijet 
mass cross section. They are the uncertainties in calorimeter calibration (cal), jet 
fragmentation (frag), underlying event (uevt), calorimeter stability over time (stab), 
relative jet energy scale as a function of pseudorapidity (rel), detector simulation 
(sim), the unsmearing procedure (unsm), the tails of the resolution function (tails), 
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absolute luminosity |T2[ of the jet 100 trigger (lum), and the relative luminosities 
of the jet 20, 50 and 70 triggers (J20, J50, and J70). The first four systematic 
uncertainties QXJ] are equivalent to a combined uncertainty in the determination of 
the dijet mass variable which decreases from 2.7% at M — 188 GeV/c 2 to 2.3% 
at M = 968 GeV/c 2 . The uncertainty in detector simulation results from a 0.5% 
uncertainty in the equality of the true dijet mass and the simulated dijet mass after 
all jet corrections are applied, independent from the first four systematic uncertainties 
mentioned above. To check that our unsmearing procedure is internally consistent, we 
applied the unsmearing procedure to a simulated dijet mass spectrum. The resulting 
Ki were in agreement with the ratios of the simulated spectrum to true spectrum for 
each mass bin. Due to limited Monte Carlo statistics, the systematic uncertainty on 
the consistency of the unsmearing procedure was 4%. The uncertainty in the dijet 
mass resolution due to non-Gaussian tails was estimated by repeating the unsmearing 
procedure with a Gaussian resolution. The systematic uncertainties on the luminosity 
for the jet 20, 50 and 70 triggers came from the statistical uncertainty in matching 
the cross section of each trigger with the next higher threshold trigger (jet 70 was 
required to match jet 100 in the first bin of the jet 100 sample, jet 50 was required to 
match jet 70, etc.) Each of the independent systematic uncertainties in Table I are 
completely correlated clS cl function of dijet mass. 

In Table II we present the fully corrected inclusive dijet mass spectrum for pp —> 
2 jets + X, where X can be anything, including additional jets. We tabulate the 
differential cross section versus the mean dijet mass in bins of width approximately 



equal to the dijet mass resolution. Figure [I] shows the fractional difference between our 



data and O(af) QCD predictions from the parton level event generator JETRAD [13 
Here the renormalization scale is \i — 0.5E™ ax , where E™ ax is the maximum jet Et 
in the generated event. In the JETRAD calculation, two partons are combined if 
they are within R sep = 1.3R, which corresponds to the separation of jets in the 
data. Predictions are shown for various choices of parton distribution functions: 



CTEQ4M and MRST @ are standard sets and CTEQ4HJ adjusts the gluon 



distribution to give a better fit to the CDF inclusive jet E T spectrum at high E T . 
Figure p] shows that the CTEQ4HJ prediction models the shape and normalization 
of our dijet data better than CTEQ4M. The CTEQ4M prediction changes by less 
than 5% when the renormalization scale is changed to /i — E™ ax , but it decreases 
between 7% and 17% for \x = 2EJp ax , and it decreases between 25% and 30% for 
fj, = 0.25E™ ax . In Fig. |2| we compare the fractional difference between our data and 
QCD with that of the D0 experiment. The D0 measurement @] and the JETRAD 
prediction obtained by D0 required that each jet be in region < 1.0. Figure |2] 
shows that our data and the D0 data are in good agreement. 

The covariance matrix for the dijet mass differential cross section is defined as 
V i:j = Sijaf(stat) + Y,f =1 G i (sys k )a j {sys k ). Here 5 {j = 1(0) for % = j(i ^ j), cr^stat) is 
the statistical uncertainty in mass bin i, and the sum is over each of the 12 systematic 
uncertainties (Ti(sysk) listed in Table I. Since the theory always predicts a smaller cross 
section than the data, the positive percent systematic uncertainty given in Table I 
was multiplied by the theoretical cross section to determine the <Ji(sysk)- From the 
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inverse of the covariance matrix, {V' 1 )^, and the difference between the data and 
the theory in each bin, Aj, we perform a x 2 comparison between the data and the 
theory. Table III presents values for % 2 = EijA^V^yAj and the corresponding 
probability for a standard \ 2 distribution with 18 degrees of freedom (14 degrees of 
freedom for the row labeled Fit). Our data is in agreement within errors with the 
QCD prediction using CTEQ4HJ parton distributions, which has an enhanced gluon 
distribution at high Et- Our data excludes CTEQ4M parton distributions, which 
have a standard gluon distribution. The x 2 comparison shows that our data cannot 
exclude with high confidence QCD predictions using MRST parton distributions, even 
though the normalization of that prediction is well beneath that of our data. This is 
because of the presence of correlated systematic uncertainties that are large compared 
with the statistical uncertainties. Such correlated uncertainties can accommodate 
certain significant deviations in both normalization and shape between the data and 
the theory with a relatively small penalty in \ 2 ■ Any theoretical prediction whose 
deviation from the data matches the shape of a correlated uncertainty will give a 
reasonable \ 2 provided that the normalization difference between the data and the 
prediction is no more than a few standard deviations. 

In conclusion, we have measured the cross section for production of two or more 
jets in the kinematic region < 2 and | cos8*\ < 2/3 as a function of dijet invariant 
mass. The data at the highest values of dijet mass are above the QCD predictions 
using standard parton distributions, similar to the excess at high Et observed in 
previous measurements of the inclusive jet Et spectrum [IJ and the total Et spec- 
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trum 0. The CDF data are described within errors by next-to-leading order QCD 
predictions using CTEQ4HJ parton distributions, and are in good agreement with a 
similar measurement from the D0 experiment. 
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Table I: Systematic uncertainties on the cross section (see text). 



Mass 
(GeV/c 2 ) 


Systematic Uncertainty on Cross Section in % 


cal 


frag 


uevt 


stab 


rel 


sim 


unsm 


tails 


lum 


J20 


J50 


J70 
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-4-1 2 
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- 7 


-4- 8 
- 7 


4 7 
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2 


4 


2 


4 


4 
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2 
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- 4 


7 
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5 


4 


5 
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2 
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8 


4 


4 8 
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5 


4 


5 
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+ 1 4 

- 9 
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4 
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7 
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8 


4 
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Table II: For each bin we list the average dijet mass, the differential cross section, 
and the statistical and total systematic uncertainty on the cross section. 



bin edge 


average M 


da/dM 


statistical 


systematic 


(GeV/c 2 ) 


(GeV/c 2 ) 


(pb/GeV/c 2 ) 


uncertainty 


uncertainty 


180 


188 


6.07xl0 2 


3.2% 


u/ 

-17 A) 


198 


207 


3.42xl0 2 


A -I Pt-f 

4.1% 


11 Q rw 


217 


228 


1.81xl0 2 


1.0% 


-4-1 Q fW 

-16 /O 


241 


252 


9.81X10 1 


-l A fW 

1.4% 


-4-1 Q rr-f 

-11% 


265 


277 


4.98X10 1 


1.8% 


-4-1 Q rw 

[J-f u/ 

-17 /O 


292 


305 


2.78X10 1 


1.1% 


-4-1 Q Crf 

-17 /0 


321 


335 


1.43X10 1 


1.4% 


~r^w u/ 

-17 /o 


353 


368 


7.41x10° 


1.9% 


~r^w u/ 
-18 /0 


388 


405 


3.83x10° 


0.9% 


-4-9 "1 rr-/ 

±11% 


427 


446 


1.89x10° 


1.2% 


-4-9 "1 rw 
Il9% 


470 


491 


9.07xl0 _1 


1.7% 


-1-99 rw 


517 


539 


4.50xl0~ 1 


2.3% 


+23 a/ 
±20% 


568 


592 


1.90X10- 1 


3.3% 


+25 ft/ 
-21 /0 


625 


652 


7.42xl0~ 2 


5.1% 


+26 ft/ 
-22 /0 


688 


716 


2.92xl0~ 2 


7.7% 


+28 ft/ 
-23 /0 


756 


784 


1.18xl0~ 2 


11% 


+30 ft/ 

-25 /0 


832 


865 


3.57xl0~ 3 


20% 


+32 W 
-26 /0 


915 


968 


9.03xl0" 4 


33% 


+34 ft/ 
-28 /° 
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Table III: \ 2 an d corresponding probability for theoretical predictions for the dijet 
mass spectrum with various choices of parton distribution functions and renormal- 
ization scales /x = DE™ ax . The row labeled Fit is the parameterization used in the 
unsmearing (see text). 



PDF 


D 


x 2 


Probability 


CTEQ4M 


0.25 


66.0 


2.2 x 10~ 7 




0.5 


48.9 


1.1 x 10" 4 




1.0 


48.1 


1.5 x 10" 4 




2.0 


52.5 


1.7 x 10~ 5 


CTEQ4HJ 


0.5 


29.8 


4.0 x 10~ 2 




1.0 


26.1 


9.8 x 10~ 2 


CTEQ3M 


0.5 


45.7 


3.3 x 10- 4 




1.0 


55.2 


1.2 x 10~ 5 


MRST 


0.5 


38.7 


3.2 x 10~ 3 




1.0 


33.5 


1.5 x lO- 2 


MRST(gt) 


0.5 


36.1 


6.9 x 10~ 3 


MRST(g|) 


0.5 


38.3 


3.5 x 10~ 3 


Fit 




20.5 


1.2 x 10- 1 
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Figure 1: The fractional difference between the measured differential cross section 
and the QCD prediction (points) as a function of dijet mass. The band is the sys- 
tematic uncertainty. The curves are the fractional difference between other QCD 
predictions, for various choices of parton distributions, and our default QCD predic- 
tion using CTEQ4M. 
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Figure 2: The difference between CDF data and QCD (solid points as shown in 
Fig. P compared to the difference between D0 data || and QCD (open points). The 
solid curves are the D0 systematic uncertainty. 
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